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Time-resolved particle image velocimetry (TRPIV) measurements are performed down- 
stream of a swept backward-facing step, with a height of 49% of the boundary-layer thick- 
ness. The results agree well qualitatively with previously reported hotwire measurements, 
though the amplitudes of the fluctuating components measured using TRPIV are higher. 
Nonetheless, the low-amplitude instabilities in the flow are fairly well resolved using TR- 
PIV. Proper orthogonal decomposition is used to study the development of the traveling 
crossflow and Tollmien-Schlichting (TS) instabilities downstream of the step and to study 
how they interact to form the large velocity spikes that ultimately lead to transition. A 
secondary mode within the traveling crossflow frequency band develops with a wavelength 
close to that of the stationary crossflow instability, so that at a certain point in the phase, 
it causes an increase in the spanwise modulation initially caused by the stationary crossflow 
mode. This increased modulation leads to an increase in the amplitude of the TS mode, 
which, itself, is highly modulated through interactions with the stationary crossflow. When 
the traveling crossflow and TS modes align in time and space, the large velocity spikes oc- 
cur. Thus, these three instabilities, which are individually of low amplitude when the spikes 
start to occur (u-m;/Ue <0.03), interact and combine to cause a large flow disturbance that 
eventually leads to transition. 


Nomenclature 
Cc chord length 
Cy pressure coefficient, C, = Dot 
) boundary layer thickness i 
Az spanwise wavelength 
Re’ unit Reynolds number 
Tu turbulence intensity 
u’,v',w' fluctuating components of velocity (u — U) 
Uae spanwise root mean square of steady disturbance velocity, U’ 
Us msit email ems temporal root mean square of u’,v’,w’ 
U,V,W time-averaged components of velocity 
U,V, W instantaneous components of velocity in the x, y and z’-directions 
Ue boundary layer edge velocity 
Use freestream velocity 
x streamwise direction 
Le direction normal to the leading edge 
Lsh number of step heights downstream of step 
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y wall-normal direction 
Zz spanwise direction (parallel to the leading edge) 


z! direction normal to x-direction 

Subscripts 

TCF band-pass filtered in traveling crossflow frequency range (80 to 200 Hz) 
TS band-pass filtered in TS frequency range (200 to 800 Hz) 


I. Introduction 


AMINAR flow control remains a promising technique for improving the fuel efficiency of aircraft in the 
Le future. In theory and in the laboratory, laminar flow control (LFC) techniques can work quite well. 
However, real-world applications result in additional difficulties that can be detrimental if not well-managed. 
For example, disruption of laminar flow can occur if the wing surface contains any small protuberances or 
surface imperfections. These imperfections can result from insect residue, rivets, bolts, steps, gaps, paint, 
and other sources. In order for LFC to effectively delay transition in an operational environment, we need 
to gain a better understanding of how surface irregularities affect transition. This is important so that we 
can specify more reasonable manufacturing tolerances. 

Stationary crossflow typically dominates the transition process in swept-wing flows. However, steps can 
act as receptivity sites for other types of disturbances. Steps can also modify the mean flow in a way that 
destabilizes those disturbances for a short region downstream of the step. 

Researchers have recently begun to study the effect of 2D steps on swept-wing transition in more detail. In 
earlier work, the results were limited to observations of the behavior of the transition front as the step height 
changed /#I2] More recently, researchers have begun to study the flow in more detail. Eppink et al Bl reported 
detailed boundary-layer measurements downstream of a backward-facing step and identified several different 
types of instabilities that were introduced by the step and modulated by the stationary crossflow vortices. 
One of the unsteady disturbances identified downstream of the step is believed to be a traveling crossflow 
instability, similar to what Tufts et all saw in their computations. Duncan et all performed hotwire 
measurements downstream of forward- and backward-facing steps to determine the effect of the steps on 
stationary crossflow instabilities. They found that the steps caused an increase in N-factor for the stationary 
crossflow, but the stationary crossflow amplitudes were very low at the step, and therefore the uncertainty 
of the N-factors was high. Tufts et all performed computations to study the interaction between stationary 
crossflow instabilities and a two-dimensional step excrescence. They found that the backward-facing step 
did not amplify the stationary crossflow modes, but they did verify the existence of a traveling instability 
in the recirculation region downstream of the step. They concluded that transition is likely caused by an 
interaction between this traveling mode and the stationary crossflow mode. Saeed et al!®! performed detailed 
hotwire measurements downstream of 2D strips placed on a swept-wing model. The strip is essentially a 
forward-facing step followed by a backward-facing step. They found that strips placed farther downstream 
generated larger amounts of unsteadiness, and concluded that the more developed stationary crossflow was 
interacting with traveling-crossflow disturbances, leading to earlier transition. 

The complicated transition process involving protuberances (such as 2D steps) apparently involves the 
interaction of multiple types of disturbances. The current research is a follow-on experiment to the hotwire 
results described by Eppink et all! and the stereo particle image velocimetry results described by Eppink and 
Yao! In the current work, we undertake time-resolved particle image velocimetry (TRPIV) measurements 
in order to gain a better understanding of the instability interactions and the breakdown mechanism leading 
to transition downstream of a backward-facing step. 


II. Previous Experimental Results 


Eppink et all and Eppink and Yad discuss the previous hotwire and PIV measurements for the 
backward-facing step (BFS) cases in detail, but we will summarize the results here as context for the follow- 
ing discussion. Hotwire measurements were performed with and without discrete roughness elements (DREs) 
on the leading edge of the model. In all cases, transition occurred far downstream and was dominated by 
stationary crossflow instabilities until the step height was increased above approximately 49% of the local 
boundary-layer thickness (6) for the baseline case. At that step height, transition moved upstream abruptly, 
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but still occurred more than 300 step heights downstream of the step. This was also well downstream of the 
separation bubble, which extended approximately 30 step heights downstream of the step. The hotwire data 
revealed a rich velocity spectrum downstream of the step, with unsteady disturbances in a broad frequency 
band (f ~ 80 to 1500 Hz). These unsteady disturbances did not occur in the baseline case and are believed 
to be directly responsible for transition in the presence of the step since the stationary crossflow amplitudes 
remained too low to cause transition via their high-frequency secondary instabilities. Using simultaneous 
measurements from a second fixed hotwire, the unsteady disturbances were categorized as three distinct 
disturbance types: a traveling-crossflow (TCF) instability (80 to 200 Hz), a Tollmien-Schlichting (TS) insta- 
bility (200 to 800 Hz), and a shear-layer instability (800 to 1500 Hz). The traveling crossflow instability was 
found to have a wavelength of approximately 50 mm, which is about 4 to 5 times larger than the predicted 
most unstable traveling crossflow mode for the no-step case. The TS instabilities were not present in the no- 
step case, and were likely amplified due to the short pressure gradient region caused by the backward-facing 
step. 

The stationary crossflow modulated the unsteady disturbances in the spanwise direction, resulting in 
regions of peaks and valleys in fluctuation amplitude with spanwise spacing corresponding to the dominant 
stationary crossflow mode wavelength. Large spikes in velocity occurred well upstream of the breakdown 
location and appeared to be the manifestation of a type of breakdown mechanism resulting from the unsteady 
disturbances. The hotwire time trace in Fig. [la] shows an example of a location which experienced a large 
number of positive velocity spikes. 
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Figure 1: Results from hotwire campaign showing measurements of the spike breakdown 
mechanism. 


Regions of positive and negative spikes occurred with spanwise spacing related to the stationary crossflow 
wavelength of 11 mm. Positive spikes occurred near the wall while the negative spikes occurred off the wall 
and offset spanwise from the positive spikes. This is illustrated in Fig. in which positive skewness indicates 
the occurrence of positive spikes, and negative skewness indicates negative spikes. Stereo PIV measurements 
revealed that these spike events were not simply local events. In fact, these events occurred at the same 
time across multiple wavelengths of the stationary crossflow instability, in a way that increased the spanwise 
modulation initially caused by the stationary crossflow instability, as illustrated in Fig. 


Ill. Experimental Setup 


The experiment was performed in the 2-Foot by 3-Foot Low Speed Boundary-Layer Channel at NASA 
Langley Research Center. The tunnel is a closed circuit facility with a 0.61-m high by 0.91-m wide by 6.1-m 
long test section. The tunnel can reach speeds up to 45 m/s (unit Reynolds number, Re’ = 2.87 x 10°/m) 


in the test section. Freestream turbulence intensity levels, Tu = Giy/3 


Flu? + v2 + w'2), were measured 


using a crosswire in an empty test section to be less than 0.06% for the entire speed range of the tunnel, and 
less than 0.05% for the test speed of 26.5 m/s. This value represents the total energy across the spectrum, 
high-pass filtered at 0.25 Hz. Thus, this tunnel can be considered a low-disturbance facility for purposes of 
conducting transition experiments!) 

The 0.0127-m thick flat plate model consists of a 0.41-m long leading-edge piece, swept at 30°, and a 
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(c) PIV data from single point (z’=-7 mm, y=0.7 mm) 


Figure 2: Instantaneous PIV snapshots at x,, ~ 120 showing spike and non-spike events. 


larger downstream piece (see Fig. [3). The model is 0.91 m wide (thus, spanning the width of the test section) 
and 2.54 m long on the longest edge. The downstream or leading edge pieces can be adjusted relative to each 
other using precision shims to create either forward-facing or backward-facing 2D steps of different heights, 
parallel with the leading edge. The leading edge piece was polished to a surface finish of 0.2 wm, and the 
larger downstream plate had a surface finish of 0.4 wm. A leading-edge contour was designed for the bottom 
side of the plate in order to make the suction-peak less severe, and therefore avoid separation, which could 
potentially cause unsteadiness in the attachment line. 

A 3D pressure body along the ceiling was designed to induce a streamwise pressure gradient, which, along 
with the sweep, causes stationary crossflow growth. A second purpose of the ceiling liner was to simulate 
infinite swept-wing flow within a mid-span measurement region of width 0.3 meters. This was achieved by 
designing the liner such that the C, contours were parallel with the leading edge within the measurement 
region. The ceiling liner was fabricated out of a hard foam using a computer-controlled milling machine. 

All measurements were performed at a freestream velocity of 26.5 m/s (Re’ = 1.69 x 10°/m), and with 
a backward-facing step height of 1.184 mm (49% 6). The current experiment utilized a single leading-edge 
roughness configuration consisting of DREs with a diameter of 4.4 mm. The DREs were applied with a 
spanwise spacing, A,, of 11 mm and were approximately 20 wm thick. The spacing of the DREs (11 mm) 
corresponds to the most amplified stationary crossflow wavelength calculated for the baseline case with no 
step. For more details of the experiment setup, refer to Eppink 

A high-speed double-pulsed Nd:YAG laser provided the laser sheet for the time-resolved particle-image 
velocimetry (TRPIV) measurements (see Fig. [3). The laser sheet was set up parallel with the leading edge 
and the backward-facing step. Two high-speed 4 megapixel cameras were used to acquire the images. One 
was placed on the outboard side of the test section at approximately 30° to the laser sheet, and the second 
camera was placed on the inboard side (in backward scattering) at an angle of approximately 45° to the laser 
sheet (Fig. |3). To achieve the desired field of view and resolution, 300 mm lenses were utilized, resulting 
in a total possible measurement area of approximately 60 mm x 30 mm. The area of interest was reduced 
to approximately 37 mm x 5 mm to obtain an acquisition rate of 4.7 kHz. This area allowed acquisition of 
approximately three wavelengths of the stationary crossflow instability in a single frame, while still acquiring 
approximately 30 points (using 75% overlap and 24x24 pixel interrogation size) inside the boundary layer. 
The cameras and laser were all mounted on the same traversing system, which allowed measurements at 
multiple locations with relative ease. An oil-based fog machine generated the seeding, which was introduced 
downstream of the test section. 
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Figure 3: Top view of PIV setup. 


IV. Data Analysis 


Data acquired at each streamwise station consisted of 10,000 images at a frame rate of 4.7 kHz. Since this 
is a stereo configuration, all three velocity components were acquired. The velocity components acquired 
from the initial processing of the images were in the direction normal to the step (#,), normal to the wall (y), 
and parallel to the step (z). The velocity components in each frame were first interpolated to remove any 
dropout points. Then the x, and z components were transformed into the components parallel and normal 
to the free stream direction (# and z’), for further analysis. 

For analysis of the two different unsteady instabilities (traveling crossflow and TS), the time series data 
at each point were filtered around the two frequency bands (80 to 200 Hz and 200 to 800 Hz) by performing 
zero-phase filtering using a 4th-order Butterworth band-pass filter. Due to the noisiness of the signal, the 
bandpass-filtered data were also smoothed in space using a 2-D Gaussian smoothing kernel with a standard 
deviation of 4 in the z-direction, and 2 in the y-direction. This smoothing was applied to aid in the detection 
and visualization of the 2D flow structures that were visible from the TRPIV data. A rectangular kernel 
was used since the flow structures are typically larger in the spanwise direction compared to the wall-normal 
direction. 

Proper orthogonal decomposition (POD) was also utilized throughout the paper to help identify flow 
structures. This method, first proposed by Lumley! is commonly used in the analysis of PIV results to 
identify the large-scale flow structures that are present in the flow. The input flow field data are deconstructed 
into a set of POD modes, which can be used to partially or fully reconstruct the flow field. The modes are 
orthogonal and are optimized such that the smallest number of modes captures the largest fraction of the 
total energy of the flow field. The modes are ranked based on their energy content, so that the first mode 
contains the most energy throughout the entire data set, and so on. 

POD is used in the current work for two main purposes. The first is to reconstruct the flow field in time 
from a small number of POD modes so that the development of the traveling crossflow and TS instabilities 
can be studied and more easily visualized (Section [VI). Without using POD reconstruction, the results are 
noisy and it is difficult to obtain useful isosurfaces, particularly when the amplitudes of the instabilities are 
low near the step. The POD reconstruction for a specific snapshot is performed by summing up the first 5 
modes multiplied by their respective coefficients for that snapshot. The second use of POD in this paper is 
to identify patterns that occur in the traveling crossflow and TS instabilities during specific flow field events 


(Section [VIp. 
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V. Comparison of TRPIV and Hotwire Data 


A comparison between the mean PIV results from the low-speed stereo PIV system and the hotwire 
results was presented by Eppink and Yao lal Overall, the PIV and hotwire results agreed well. In this section, 
we will focus on a comparison of the fluctuating velocity component. Velocity spectra acquired at several 
different streamwise locations downstream of the step at y=1 mm are shown in Fig.|4}| These spectra were 
averaged across the span in order to get a sense of the overall unsteady flowfield. Results from the hotwire 
data are shown in Fig. while results from similar streamwise locations are shown in Fig. for the 
TRPIV measurements. Figure also includes a single hotwire spectra (from x2,,=69) so that a direct 
comparison can be made. Note that the y-scales are different between the two plots. The TRPIV data 
have a much higher noise floor than the hotwire data (about 5 to 6 orders of magnitude), which is to be 
expected. However, we can still clearly discern two broad peaks between 80 to 200 Hz and 200 to 800 Hz in 
the TRPIV spectra corresponding to the traveling crossflow and TS frequency bands. These two frequency 
ranges are denoted by the vertical black lines in both figures. In general, the TRPIV spectra are larger in 
amplitude compared to the hotwire spectra from similar locations. This is especially true for the locations 
closer to the step, where the amplitudes are lower, such as x5, ~ 20. Notice also that there is not very much 
energy present in the higher frequency band (800 to 1500 Hz) in the TRPIV spectra, which was associated 
with the shear-layer instability from the hotwire data. Despite the discrepancy in amplitudes and the lack 
of resolution of the higher-frequency modes, it is quite encouraging that we are even able to measure these 
fluctuations at all using TRPIV, since they are of fairly low amplitude. 

One question that arises when performing TRPIV is the ability of the particles to respond to the high 
frequency perturbations in the flow. Using the approach of Mei #4) we can estimate a low-pass frequency 
cut-off of the particles, assuming a certain particle diameter and density ratio. If we assume the particle 
density is approximately that of water, and the particle diameter ranges from 0.5 to 2 ym, the estimated 
low-pass frequency cutoff ranges from 14 to 227 kHz. Thus, the particles in this flow should be able to 
respond to frequencies in the 800 to 1500 Hz range. 
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Figure 4: Comparison of spanwise-averaged velocity spectra from hotwire and TRPIV results 
at several streamwise locations downstream of the step. 


The TRPIV and hotwire results are further compared in Fig.[5| in which color contours of u/,,,,, amplitudes 


are plotted within the traveling crossflow and TS frequency bands across the full y-z plane. Unfortunately, 
we do not have full planes of hotwire data and TRPIV data for the same DRE case at the same streamwise 
location, but we can still qualitatively compare results from different streamwise locations. Fig. |5}| shows 
TRPIV data at 2,,=73 and hotwire data at x,,=100. The amplitudes for the TRPIV data are higher than 
that of the hotwire data, even though the TRPIV data were acquired upstream of the hotwire data, where 
amplitudes should be lower. This agrees with the general trend from the spectra. However, we still observe 
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the same modulation patterns due to the interaction of the stationary crossflow with the different unsteady 
instabilities. The w/.,, contours for the traveling crossflow frequency band exhibit a modulation pattern 
with a wavelength of half of the primary stationary crossflow mode wavelength of 11 mm, with some peaks 
in amplitude occurring near the wall, and the adjacent peaks occurring shehely off of the wall. This pattern 
is similar between the TRPIV and hotwire data (Figs. [Bal and [Sch . The u/,,,, contours for the TS frequency 
band in both cases (Figs. [5b] and [5dp are strongly modulated eal a wavelength of 11 mm, equal to that of 
the primary stationary crossflow mode. Overall, the qualitative agreement between the TRPIV and hotwire 
data is encouraging. Despite the discrepancy in amplitude of the fluctuating component, the results can still 
provide valuable insight into the instability interactions that occur downstream of the step. 
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Figure 5: Comparison of fluctuating velocity RMS for traveling crossflow and TS frequency 
bands. 


VI. Development of Unsteady Instabilities 


In this section, the downstream development of the traveling crossflow and TS disturbances is examined. 
As mentioned in the previous section, the higher-frequency band associated with the shear-layer instability 
was not well-resolved with the TRPIV system due to the higher noise level. Additionally, this instability is 
of lower amplitude than the TCF and TS instabilities, and therefore is not believed to be as important in the 
transition process downstream of the step. Due to these reasons, analysis of this higher-frequency instability 
is not included in this paper, and the discussion instead focuses on the two lower frequency instabilities. 

Data were acquired at several stations downstream of the step, starting as close to the step as physically 
possible, given the finite thickness of the laser sheet and the inevitable reflections that occur near the step. 
This location is estimated to be approximately 1 mm downstream of the step, though it is not possible to 
obtain an exact location since we do not know precisely where the center of the laser sheet is located. The 
fluctuation amplitudes at this station were so low compared to the noise level that there was nothing of 
interest to report. Therefore, the measurements reported in this paper begin at x,,=9 and go downstream 
to 2s,=73 in approximately 13 mm increments. 

POD was performed individually on the traveling crossflow and TS band-pass-filtered data to examine 
the dominant flow structures of each instability. The POD energy spectra for the first 50 POD modes at all 
locations are shown in Fig. [6] for the traveling crossflow and TS data. The eigenvalues of the POD modes 
correspond to the kinetic energy content of the respective POD mode. Thus, the cumulative modal energy 
corresponds to the cumulative sum of the eigenvalues, normalized by the sum of all the eigenvalues. The 
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plots indicate that at most stations (excluding the most upstream station), the first mode makes up at least 
40% of the total modal energy. Moreover, the first 5 POD modes account for more than 90% of the total 
modal energy. The lower values of cumulative energy for the most upstream station («,,=9) are believed 
to be due to the extremely low signal-to-noise ratio at this station, since the amplitudes of the disturbances 
were barely discernible. 
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Figure 6: Cumulative POD eigenvalues (energy) for the first 50 POD modes at each streamwise 
location. 


The first POD modes at each streamwise station for the traveling crossflow band-pass-filtered data are 
shown in Fig. |7|for all three components of velocity. The first POD modes for the three most upstream 
stations consist of a structure with a large spanwise wavelength of approximately 40 mm. In the hotwire 
data presented by Eppink® the spanwise wavelength of the traveling crossflow was found to be approximately 
50 mm. The TRPIV measurement is less than the 50 mm wavelength found from the hotwire data, but 
still in good qualitative agreement with the findings. The stations downstream of reattachment, starting at 
Lspn=41, begin to exhibit some modulation of that primary large wavelength structure. Eventually, at the 
last two stations (x,,=62 and 73), the first POD mode exhibits a strong 11 mm wavelength, and the large 
wavelength structure is no longer present. 

The first POD modes at each streamwise station for the TS band-pass filtered data are shown in Fig. 
The first mode at the most upstream station (%sn=9) for the TS band-pass-filtered data is mostly dominated 
by noise near the left edge of the figure. However, all of the stations downstream show a clear modulation of 
the TS instability by the primary stationary crossflow mode. The first POD mode for this instability does 
not change drastically downstream, but the modulation becomes stronger and the peaks more localized. 
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Figure 7: First POD modes at each streamwise station for traveling 
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Figure 8: First POD modes at each streamwise station for TS frequency band. 


Figure. [9] shows band-pass-filtered results at three sequential times for the traveling-crossflow frequency 
band at 2,,=19. The right column of this figure includes the reconstructed data at the same times using the 
first five POD modes. While the large traveling crossflow structure is visible in the bandpass-filtered data 
(left column), it is very noisy. The POD reconstruction (right column) provides a better visualization of 
the behavior of this instability. This instability travels from left to right primarily in the spanwise direction 
(in-plane), which is in agreement with the hotwire measurements. 

Figure shows the same results but for the TS frequency band. Once again, the TS structures are 
somewhat visible in the band-pass-filtered data, but the patterns become much more clear with the POD 
reconstruction. This mode, in contrast to the traveling crossflow mode, is traveling primarily out of plane. 
Therefore, there is not much in-plane motion visible, particularly at the most upstream stations, before it 
becomes strongly modulated by the stationary crossflow vortices. 
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Figure 9: Band-pass filtered streamwise velocity in the traveling crossflow frequency range 
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Figure 10: Band-pass filtered streamwise velocity in the Tollmien-Schlichting frequency range 
(200 to 800 Hz) at x,,=19 for three consecutive times. Figures in the left column are the 
bandpass-filtered data, the figures in the right column are the POD reconstruction of the 


bandpass-filtered data using the first 5 modes. 
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We can further visualize the behavior of these two instabilities by expanding the results in three dimen- 
sions, using time as the x-axis, and plotting isosurfaces of the band-pass-filtered streamwise velocity. POD 
reconstruction using the first 5 modes is performed for 100 samples (constituting a total of 0.02 s). The 
POD-reconstructed data are used here because the band-pass-filtered data are too noisy to obtain useful 
isosurfaces, particularly for the most upstream stations. The results are shown in Fig. [11] for the first three 
stations, and Fig. for the last four stations. Note that these plots can be somewhat misleading, par- 
ticularly for the traveling crossflow instability, since this instability is traveling primarily in plane. Thus, 
these are not the actual three-dimensional spatial structures that exist, this is simply a convenient way of 
visualizing the action of this instability over time in the 2D plane. 

It should be noted that we cannot draw firm conclusions about the physical nature of the disturbances 
solely from these data since the data were only acquired in a single plane. Thus, it is not clear which direction 
the instabilities are traveling. However, we know from the hotwire results that, at least at the most upstream 
stations before the instabilities become highly modulated by the stationary crossflow, the traveling crossflow 
mode is traveling primarily in the spanwise plane, while the TS mode is traveling close to the chordwise 
direction. 
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Figure 11: Traveling crossflow and TS streamwise velocity isosurfaces reconstructed using the 
first 5 POD modes and expanded in time. First three stations shown from r,,=9 to 30. 


The traveling crossflow mode (left column) appears as a spanwise traveling disturbance as observed from 
Fig. [9] at the first two stations. Starting at x,,=30 (Fig. (11), some modulation of the structures starts 
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to become visible at the 11 mm primary stationary crossflow mode wavelength. This modulation becomes 
stronger and stronger downstream, until 7,,=62 (Fig.[12), where the large spanwise traveling structure from 
the upstream stations is no longer readily visible. The TS mode (right column) appears as a 2D, streamwise- 
traveling mode for the first three to four stations. Starting at x,,=19, some minor spanwise modulation 
of the structure starts to become visible. This modulation becomes stronger downstream. By x5,=52, the 
mode has become highly three-dimensional. In fact, it looks like there may be two oblique modes present, 
though it is not possible to tell from these results alone. Wasserman and Kloker2l studied the interaction 
of TS and stationary crossflow instabilities in a flow with a pressure-gradient changeover. They found that 
the interaction resulted in the growth of a pair of oblique modes, which they referred to as low-frequency 
secondary modes. They also found that these modes experienced localized amplification in the spanwise 
plane due to the deformed primary state. We may be seeing a similar phenomenon in this flow. 
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Figure 12: Traveling crossflow and TS streamwise velocity isosurfaces reconstructed using the 
first 5 POD modes and expanded in time. Last four stations shown from 2,,=41 to 73. 
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VII. Velocity Spikes 


The aforementioned large amplitude velocity spikes first become apparent at the x,,=52 station. The data 
at this station are therefore analyzed further to attempt to gain a better understanding of this mechanism. 
The time-averaged streamwise velocity field is shown in Fig.|13/along with the skewness of the time-resolved 
data. The skewness, similar to the hotwire data (Fig. rip divorces that positive spikes (indicated by 
positive skewness) occur mainly near the wall, while negative spikes occur off the wall and slightly offset 
in i} from the locations of positive spikes. The time series from a single point (z=-3 mm, y=0.5 mm in 
Fig. bas is shown in Fig. [14] [14] This point is located in a region where positive velocity spikes occur. 
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(a) Time-averaged streamwise velocity at 75,=52. (b) Skewness of streamwise velocity at 1;,=52. 


Figure 13: Mean flow and skewness at 2,;,=52. 
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Figure 14: Time trace at 7,,=52, y=0.5 mm, z=-3 mm. 


A single snapshot during one of the large velocity spikes is shown in Fig. which includes the stream- 
wise velocity along with the band-pass-filtered streamwise velocity for both the traveling crossflow and TS 
frequency bands. At this streamwise location, the stationary crossflow has caused a strong modulation of the 
traveling crossflow, which is evident in Fig. During the spike event, the traveling crossflow has reached 
a particular point in the phase such that it is positive near the wall and negative away from the wall. These 
positive and negative regions occur in certain spanwise locations such that the overall effect is an enhance- 
ment of the modulation caused by the stationary crossflow. The stationary crossflow also modulates the TS 
instability, but in a different manner. At this streamwise location, the initially 2D TS wavefront has become 
distorted such that the negative and positive parts of the phase occur at the same time in this spanwise 
plane. The motion (which is illustrated by Eppink et al 23) using phase-locked hotwire results), is such that 
the positive peak occurs off the wall, then moves toward the wall and to the left. Then the negative peak 
does the same. 

Note that there are regions where the positive parts of the traveling crossflow and TS instabilities occur 
simultaneously (such as near the wall at z ~ -3). The same is true of the negative parts (such as z & 0, 
y =1). These correspond to the regions where the positive and negative spikes occur, respectively. In other 
words, when the phases of the traveling crossflow and TS instabilities align in time and space, the large 
spikes in velocity occur. However, when the traveling crossflow reaches the opposite part of its cycle (i-e., 
180° out of phase with Fig.|15b), such that it counteracts the spanwise modulation caused by the stationary 
crossflow, the spikes do not occur, even when the TS and traveling crossflow instabilities align in space. 
This is why we only observe positive and negative spikes in certain spatial areas. To further confirm these 
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observations, we perform POD analysis. First, we will look at the POD modes for the entire set of 10,000 
snapshots. Second, we will perform the POD analysis on a subset of these snapshots that correspond to the 


times when the large velocity spikes occur. 
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Figure 15: Streamwise velocity, TCF and TS band-pass filtered data during a spike event at 
Lsgn—O2. 


The POD analysis is performed individually on the two band-pass-filtered data sets. The first five 
modes from the POD analysis on the full data set are shown in Fig. At this streamwise location, the 
instabilities are beginning to undergo a strong spanwise modulation due to the action of the stationary 
crossflow vortices. All five POD modes of the TS band-pass-filtered data (right column) exhibit a clear 
11 mm spanwise wavelength. The second mode (Fig. looks similar to the TS snapshot during the 
spike event (Fig. |15c). The first two modes of the traveling crossflow band-pass-filtered data (Figs. 
and show some remnant of the large wavelength spanwise traveling mode that was evident at the 
upstream stations (Fig. [9). The shift in phase between mode one and two indicates that these modes form 
a pair, physically indicating the presence of a traveling mode in the spanwise plane. The third through 
fifth modes (Fig. [16e]g,i), however, exhibit a clear spanwise wavelength of the dominant stationary crossflow 
mode (11 mm). The third mode is similar to the TCF snapshot during the spike event (Fig. |15b). This 
mode, when superposed onto the mean flow with a positive coefficient, results in an increase in the spanwise 
modulation of the flow initially caused by the stationary crossflow. This effect is illustrated in Fig. which 
shows the time-averaged velocity at y=1 across the span compared to the time-averaged velocity with the 
addition of the traveling crossflow POD mode 3 (Fig. with a large positive and negative coefficient. As 
this figure illustrates, when this mode is superposed with a positive coefficient, there is an increase in the 
spanwise modulation compared to the mean flow. However, the opposite is true when a negative coefficient 
(representing a 180° phase shift) is applied. 
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Figure 16: First five POD modes for the traveling crossflow and TS bandpass-filtered full data 
sets, ©.,=52. 
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Figure 17: Comparison of spanwise modulation of mean flow with and without TCF POD 
mode 3. 


We explore the traveling crossflow modes further to try to determine the physical nature of the different 
sets of POD modes. Since longer wavelength structures are apparent in the first two POD modes, and the 11 
mm wavelength structures are more dominant in modes 3 through 5, it appears that these two sets of modes 
may constitute different physical structures. A 0.085 second long stretch of data is reconstructed using just 
the first two POD modes (Fig. [18ah, and then separately the third through fifth POD modes (Fig. [18b). As 
expected, the reconstruction using the first two modes shows a slightly modulated spanwise-traveling wave 
with a long wavelength, similar to that seen in the POD mode reconstructions farther upstream (Fig.[11). The 
reconstruction using modes 3 through 5 illustrates a very different structure that exists in this same frequency 
band. This structure may also be a spanwise-traveling wave, but with a smaller spanwise wavelength, A, = 
8.5 to 20 mm, calculated based on the phase speed in the spanwise plane and the frequency band of 80 to 
200 Hz. During this particular span of time, the mode primarily appears to be traveling in the outboard 
(positive z) direction, though there are times where it changes direction and appears to be traveling in the 
opposite direction, such as t=0 to 0.01 s. It is also possible that this mode is not a spanwise traveling 
crossflow instability, but instead an oblique wave, or a pair of oblique waves similar to that seen in the TS 
frequency band starting at 7.,=52 (see Fig. (12). Perhaps this is another set of oblique TS modes resulting 
from the TS-stationary crossflow interaction, but at a lower frequency. Unfortunately, we cannot draw any 
definitive conclusions about what type of instability this is based on the available data. 

Next, the POD analysis is applied to the subset of snapshots corresponding to the times when the large 
velocity spikes occur. These samples are identified by using the time series shown in Fig. A threshold 
of u’=1.7 m/s is chosen to be representative of a spike event (which is twice the standard deviation of wu’ (t) 
at this location), and the peaks above this threshold are identified from the time series. This results in 108 
spike samples over which to perform POD. Note that we have also performed this analysis using a time series 
from a different location (such as a location where negative spikes occur), and the results are similar since 
the spikes occur nearly simultaneously across the whole plane of data. The first POD modes for the traveling 
crossflow and TS frequency bands from this analysis are shown in Fig. These two modes also appear in 
the POD modes from the full set of data (Fig. (16), but they are not the first modes. The traveling-crossflow 
first mode from the POD analysis of the spike samples (Fig. corresponds to the third mode from the 
POD analysis of the full data set (Fig. [16c). The TS first mode from the POD analysis of the spike samples 
(Fig. [19b) most closely resembles the second mode from the POD analysis of the full data set (Fig. [16d). 
The POD energy spectra for the first 50 POD modes for both the traveling crossflow and TS results are 
shown in Fig. The first mode for the TCF data contains over 50% of the total energy, while the first 
mode for the TS data contains over 70% of the total energy. Thus, these modes are clearly very dominant 
during the spike events. 
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Figure 18: Reconstruction of TCF frequency band using first two POD modes and modes 3 
through 5 separately, 7,,=52. 
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Figure 19: First POD mode from spike subset of data. 


To illustrate the effect of the spike events on the entire flow field, we also show the first POD mode from 
the spike subset of snapshots taken from the unfiltered data set. This first POD mode is shown in Fig. 
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and essentially looks like a combination of the TCF and TS first modes shown in Fig. Figure 21b] shows 
a superposition of the time-averaged streamwise velocity (shown in Fig. with the POD mode applied 
with a coefficient of 50, showing the effect of this POD mode on the mean flow. The enhanced modulation 
of this flow field is evident when compared to the mean flow in Fig. [13a] and is similar to the effect seen for 
the single snapshot during the spike event shown in Fig. 
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Figure 20: Cumulative modal energy of first 50 POD modes from spike subset of data. 
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Figure 21: First POD mode from spike subset of unfiltered data and its effect on the mean 
flow. 


Histograms of the POD coefficients for the first modes from the spike subset of data are shown in 
Fig. The coefficients are mostly positive for both modes, indicating that the spikes primarily occur 
during the positive part of the phase (i-e., when they are both positive near the wall, as illustrated in Fig. [19p. 
Additionally, a scatter plot of the coefficients for the traveling crossflow and TS first modes (Fig.[23) indicates 
a direct relationship between the amplitude and phase (i.e, sign) of the two modes during the spike events. 
In context of the full data set, these two modes do not occur more often with a positive phase rather than a 
negative phase. This can be illustrated by examining the histograms of the POD coefficients of the modes 
from the full data set that most closely resemble the first modes from the spike subset. Fig [24] shows the 
histograms of the coefficients for the third mode from the traveling-crossflow-filtered data (Fig. |16e), and 
the second mode from the TS-filtered data (Fig. [16d). In contrast to the coefficients from the spike subset of 
data (Fig. (22h, these histograms are centered around zero and are not skewed either positively or negatively. 
This indicates that the positive phase of these modes is correlated with the spike events, but in general the 
positive and negative phases occur equally. 
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Figure 22: First POD mode coefficients from spike subset of data. 
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Figure 23: Scatter plot of first mode coefficients from spike subset of data. 
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Figure 24: POD mode coefficients from full data set. 
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Scatter plots of these POD coefficients are shown in Fig. [25]in order to examine the relationship between 
these TCF and TS modes. We plot the first POD mode for the TS-filtered data in addition to the second 
POD mode because these two modes make up the majority of the energy for this instability. In both figures, 
the x-axis is the coefficient of the third POD mode of the traveling-crossflow filtered data. This value is 
plotted against the first mode coefficient for the TS filtered data in Fig. and against the second mode 
coefficient in Fig. |25b| These plots show an interesting relationship between these two sets of coefficients. 
Specifically, as the coefficient of the traveling-crossflow mode increases, the scatter of both of the TS POD 
modes increases (i.e., they reach larger absolute values). To further illustrate this relationship, we divided 
the TCF mode coefficients into 20 bins and calculated the rms amplitude of the TS POD mode coefficients 
within each bin. These values, multiplied by +2 to illustrate the peak amplitudes, are shown as the red 
lines in these figures. These lines clearly illustrate the effect of the sign and amplitude of the TCF 3rd POD 
mode coefficient on the fluctuation amplitude of the two TS POD mode coefficients. This behavior indicates 
that the amplitude of the TS fluctuations depends on the phase and amplitude of this particular traveling 
crossflow mode. As was mentioned previously, this traveling crossflow mode causes the modulation of the 
streamwise velocity to increase during the positive part of the phase. Thus, the enhanced modulation of the 
flow may be causing an increase in the amplitude of the TS mode. 
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Figure 25: Scatter plot showing relationship between TCF mode 8 coefficients, and TS mode 
1 and mode 2 coefficients. 


The first POD traveling-crossflow modes at the two downstream stations (#,,=62 and 73 in Fig. 
look similar to the third traveling-crossflow POD mode at 2,,=52 (Fig. |16e). This indicates that this 
mode becomes more dominant at the downstream locations. Additionally, the large velocity spikes at these 
stations occur more frequently and are becoming larger in amplitude. By the last station, the flow is starting 
to become intermittent. The high frequencies typically occur during the spike events, indicating that the 
breakdown is occurring intermittently as a result of the large flow modulation that occurs during the large 
velocity spikes. Figure[26a|shows a time trace from a location near the wall at 7,,=73, zoomed in to show the 
high frequencies occurring during the large velocity spike. The instantaneous streamwise velocity field in the 
spanwise plane during the spike at t=0.231 seconds is shown in Fig. [26b] This flow field, compared to that 
shown in Fig. [15a at Lsn=52, is fairly chaotic and reveals the existence of large-amplitude small-wavelength 
structures, indicating the beginning of breakdown. 
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Figure 26: Time trace and streamwise velocity field showing spike event, 7,,=73. 


VIII. Summary and Conclusions 


Time-resolved PIV measurements of the flow instabilities downstream of a backward-facing step in a swept 
wing flow are presented. The TRPIV measurements, though much noisier than the previously reported 
hotwire measurements, were still able to resolve the low-amplitude instabilities in the flow. The TRPIV 
measurements agree well qualitatively with the hotwire results. The spectra show the same broadband 
peaks around the traveling crossflow and TS frequency band, and the contours of RMS fluctuating velocity 
reveal the same modulation patterns that occur due to interactions with the stationary crossflow vortices. 
The fluctuating amplitudes measured using TRPIV were somewhat higher than the corresponding hotwire 
amplitudes, which is likely a consequence of the additional noise inherent in the technique. However, it is 
encouraging that TRPIV was able to resolve the instabilities at such low amplitudes. The results provide 
valuable insight into the development and interactions of the instabilities that cannot be obtained using a 
point measurement technique such as hotwire anemometry. 

In the previously-reported hotwire measurements, large velocity spikes were observed to occur prior to 
breakdown downstream of the backward-facing step. These velocity spikes grew larger in amplitude down- 
stream and eventually experienced high-frequency oscillations, indicating the beginning of breakdown. The 
TRPIV results allow us to study the downstream development of the traveling crossflow and TS instabilities 
and to see how they interact to produce this particular path to breakdown. 

The TS instability becomes highly modulated due to interactions with the stationary crossflow instability. 
This was known from the hotwire data, but the TRPIV results also reveal more details about the highly 
three-dimensional nature of this instability downstream of reattachment, as well as the possibility of the 
existence of a pair of oblique modes starting at x,,=52. 

The traveling crossflow mode begins as a large wavelength (A, ~ 40 mm) spanwise traveling mode, 
which was also observed from the hotwire measurements. POD analysis indicates the existence of another 
mode within this same frequency range beginning downstream of reattachment. At a certain point in the 
phase, this mode enhances the streamwise velocity modulation initially caused by the stationary crossflow 
instability. This occurs because the spanwise wavelength of this traveling crossflow mode is similar to that 
of the primary stationary crossflow mode. It is at this point that the large amplitude velocity spikes tend to 
occur, but only when the TS mode is also at a similar point in its phase, such that the positive and negative 
parts of the mode mostly align with the TCF mode, further enhancing the modulation of the flow. However, 
the velocity spikes are not simply a superposition of the two modes. This is evident from the fact that the 
spikes do not occur during the opposite part of the phases (180° phase shift) of the TCF and TS modes. 
The POD analysis also indicates that the enhanced modulation of the flow caused by this TCF mode results 
in an increase in the fluctuation amplitude of the TS mode. Thus, the large velocity spikes appear to be a 
result of multiple interactions. First, the stationary crossflow instability causes modulations of the TCF and 
TS modes. The TCF mode, at the right point in the phase, causes an instantaneous increase in the spanwise 
flow modulation, which results in an increase in amplitude of the TS mode. The TS mode and the TCF mode 
then combine, at the right point in the phase of both disturbances, to cause a large spike in velocity. These 
large spikes cause a very large modulation of the flow, which ultimately results in intermittent breakdown. 

These results highlight the mechanism through which the instabilities in a low-amplitude broadband 
disturbance field in a three-dimensional flow can interact to cause large amplitude disturbances, and ulti- 
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mately, breakdown and transition. The fluctuating instabilities, whether initially 2D or 3D, become highly 
modulated through interactions with the stationary crossflow vortices, resulting in points in the phase where 
they increase the spanwise shear of the flow at certain locations. This further destabilizes the instabilities at 
these locations, leading to even larger modulations of the flow, eventually resulting in the large amplitude 
velocity spikes when the multiple types of instabilities are superposed at these locations. In this scenario, 
the stationary crossflow instabilities, though not strongly impacted by the backward-facing step, play an 
important role in transition through the modulation of the mean flow and the unsteady instabilities. 
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